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ABSTRACT 

Aims. We report the X-ray properties of young (~3 Myr) bona fide brown dwarfs of the Taurus Molecular Cloud (TMC). 
Metliods. The XMM-Newton Extended Sun'ey of the TMC (XEST) is a large program designed to systematically investigate the X-ray 
properties of young stellar/substellar objects in the TMC. In particular, the area surveyed by 15 XMM-Newton pointings (of which 
three are archival observations), supplemented with one archival Chandra observation, allows us to study 17 brown dwarfs with M 
spectral types. 

Results. Half of this sample (9 out of 17 brown dwarfs) is detected; 7 brown dwarfs are detected here for the first time in X-rays. 
We observed a flare from one brown dwarf. We confirm several previous findings on brown dwarf X-ray activity: a log-log relation 
between X-ray and bolometric luminosity for stars (with L,<10Lo) and brown dwarfs detected in X-rays, which is consistent with a 
mean X-ray fractional luminosity < log(Lx/L,) > = -3.5±0.4; for the XEST brown dwarfs, the median of log(Lx/L,) (including upper 
limits) is -4.0; a shallow log-log relation between X-ray fractional luminosity and mass; a log-log relation between X-ray fractional 
luminosity and effective temperature; a log-log relation between X-ray surface flux and effective temperature. We find no significant 
log-log correlation between the X-ray fractional luminosity and £'W(Hc!'). Accreting and nonaccreting brown dwarfs have a similar 
X-ray fractional luminosity. The median X-ray fractional luminosity of nonaccreting brown dwarfs is about 4 times lower than the 
mean saturation value for rapidly rotating low-mass field stars. Our TMC brown dwarfs have higher X-ray fractional luminosity than 
brown dwarfs in the Chandra Orion Ultradeep Project. 

Conclusions. The X-ray fractional luminosity declines from low-mass stars to M-type brown dwarfs, and as a sample, the brown 
dwarfs are less efficient X-ray emitters than low-mass stars. We thus conclude that while the brown dwarf atmospheres observed here 
are mostly warm enough to sustain coronal activity, a trend is seen that may indicate its gradual decline due to the drop in photospheric 
ionization degree. 

Key words. Stars: low-mass, brown dwarfs - X-rays: stars - ISM: individual objects: the Taurus Molecular Cloud 



1. Introduction 

Pre-main-sequence low-mass stars, i.e. T Tauri stars, show a high 
level of X-ray emission, which is generally attributed to an active 
corona, an enhanced version of the magnetic corona on t he con - 
temporary Sun (see for a review Feiaelson & Montmerle j 1 999h . 
The dynamo mechanism producing the magnetic field in these 
fully convective stars is still discussed (Preibisch et al. 2005a). 
X-ray observations of brown dwarfs (BDs) allow investigation 
of the magnetic activity in the substellar regime, where masses 
are lower than about 0.075 Mq. 

The first X-ray detectio n of a BD was reported by 
iNeuhauser & Co meronI Jj998|), who identified in the Cha I dark 
cloud an X-ray source detected by ROSAT - in early 1991, be- 
fore the near-infr ared observations of the first BDs were reported 
(Nakaiir na et alj fr995: Rebolo et al. 1995) - as the counterpart 
of a young BD. ROSAT detected only a few BDs or very low- 
mass stars, all of them young and located in star-forming regions 



jNeuhauser et al JFTO IcWron et al Jl200(l iMokler & Stelzed 
l20(m 

The new generation of X-ray satellites, XMM-Newton and 
Chandra, allows now detection of more BDs thanks to 
their increase of sensitivity. The number of young or 
intermediate-age BDs detected in X-rays 
2001; Preibisch & Zinneck e3l200 L l2002t 



Briaas&Pve 2004^ .Stelzer et al. 



(e.e.. Ilmanishi et alJ 


iTsuboi et alJ 


20031 


lOzawa et alJ 


20051 
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Stelzer et al. 2006) exceeds largely the one of (older) field BDs, 
whi ch have b een detected so far mainly during X-ray flares 
(Rutledge et aL.200Cl: .Stelzer. . 2004). Recently, t he X-ray prop- 
erties of 33 young BDs dSlesnick et alJ I200I in the Orion 
Nebula Cluster (ONC) located at ~450pc were studied by 
Preibisch et al, 1 2005 h) in the Chandra Orion Ultradeep Project 
(COUP). 

We investigate in this paper the X-ray properties of a sam- 
ple of 17 BDs in the nearby (140pc) Taurus Molecular Cloud 
(TMC), surv eyed in the XMM-N ewton Extended Survey of the 
TMC (XEST lGiidel et alJ2006aft . First, we define in Clthe TMC 
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Fig. 1 . The XMM-Newton Extended Survey of the Taurus Molecular Cloud (XEST). Contours show the '^CO emission ('Da me et"al] 
Il987h of the Taurus Molecular Cloud (TMC), overlaid o n a visual extinct ion map (linear colour scale, with black colour correspond- 
ing to A\ ~ 6 mag; based on DSS I, 6'-resolution map. lDobashi et all20 05). The 27 XMM-Newton fields of view of the XEST are 
plotted with continuous/dashed circles (of which seven are archival observations and one is a separate program on T Tau; dashed 
circles). Note the outlying XMM-Newton fields of view around SU Aur (NE corner) and L1551 (S). The two squares near labels 3 
and 8 show the archival Chandra fields of view used in this work ( see Table |AJJi. The continuous line shows the survey of brown 
dwarfs (BDs) performed with the Canada-France-Hawaii telescope ("Gu ieu et aL2 006). This region hosts 42 young BDs of the TMC 
(Briceno et al. 1998, 2002; Luhman 2000, 2004, 2006; Martin et al. 2001t lGuieu e t al. 2006). Black dots indicate the 25 BDs of this 
region which were not surveyed in X-rays with the XEST. White dots show the 8 BDs not detected in X-rays. Crosses show the 9 
BDs detected in X-rays (T abled with labeled numbers referring to B D names. Only 2 BD s of TMC were previously detected in 
X-ray s by ROSAT: MHO 4 jCarkner et allT996tlNeuhauser et al.ll999h and CFHT-BD-Tau 4 jBriceno et alfT999tlMokler & Stelzeil 



BD sample, review previous X-ray results on it, and report the 
BDs detected in the XEST. Analysis of the X-ray variability of 
the TMC BDs is presented in ^ and X-ray spectral proper- 
ties are determined in ^ In ^ we compare the X-ray lumi- 
nosities of the BDs and low-mass stars observed in the XEST. 
Comparison of the X-ray fractional luminosities of BDs in the 



XEST and COUP is made in ^ We discuss the origin of the 
BD X-ray emission in the broader context of the X-ray emission 
of cool stars in ^ Finally, we summarize our results in ^ 
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Table 1. The sample of the 17 young BDs surveyed in the XEST. Col. (2) gives the reference of the discovery paper Col. (3) 
give the 2MASS counterparts. Col. (7) gives the references for the spectral type (Col. 4), the optical extinction (when not avail- 
able derived from Aj using Rieke & Lebofskv 1985, Col. 5), and the reference of the equivalent width of Ha (EW(Ha), given in 
Col. 6) if it is different. Negligible optical extinction in Col. (5) are indicated by 0. The effective temperature in Col. (8) has been 
computed from t he spe ctral type using the temperature scale T^ff - 3841.94 - 141.17 x S pTyp, which is valid for young stars with 
M spectral type dGuieu et al. 2006). The visual extinction in Col. (5) is taken when available from the literature; for 2MASS J0421 
and 2MASS J0422, we averaged the visual extinction computed from J - H, H - K^, and the spectral type, using a dwarf se- 
quence (compiled from the literature; e.g., Leagett et al. 1998). T he bolom etric luminosity of the substellar photosphere in Col. (9) 
has been computed from /, /-band magnitudes and A y (see iGuieu et alj|2006l): for 2MASS J 0414 and MHO 4 (without /-band 
magnitude available) the reference for the luminosity is'Luhman'(^2004') and'Briceno et alJ ( 12002) . respectively. Col. (10) indicates 
accreting sources based mainly on ^^(Ha) (see ^6.21 . The last column indicates detection in X-rays (thi s work; see Table 12} . 
References: B98=Briceno et al. ( 1998); B02=Briceno et al. (2002); G06 ^Guieu et al.1 ( 12003) : L04 jLuhmar3 ( l2004 : L06= lLuhrnanl 
aOOQ : MOl ^Martin et al..c2001.) : MOS ^MuzeroUe et al. (2005a). 



Name 


Ref. 


2MASS 


SpTyp 


Av 


EW(Ra) 


Ref. 




L, 


Acc. 


X 










mag 


A 




K 


Lo 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


2MASS J0414 


L04 


J04141 188+281 1535 


M6.25 


1.1 


250.0 


M05 


2960 


0.015 


y 


n 


KPN0-Tau2 


B02 


J04185115-H2814332 


M6.75 


0.4 


8.4 


G06 


2889 


0.007 


n 


n 


2MASSJ0421 


L06 


J04215450-I-2652315 


M8.5 


3.0 




L06 


2642 


0.003 




n 


2MASS J0422 


L06 


J0422 1332+ 1934392 


M8 


1.0 




L06 


2713 


0.017 




y 


KPNO-Tau4 


B02 


J04272799+26 12052 


M9.5 


2.5 


158.1 


G06 


2501 


0.004 


y 


n 


KPNO-Tau5 


B02 


J04294568+2630468 


M7.5 





30.0 


B02 


2783 


0.023 


n 


y 


KPN0-Tau6 


B02 


J04300724+2608207 


M9 


0.9 


207.9 


G06 


2571 


0.003 


y 


n 


KPN0-Tau7 


B02 


J043057 18+2556394 


M8.25 





300.0 


B02 


2677 


0.004 


y 


n 


MHO 4 


B98 


J043 12405 +18002 15 


M7 


0.5 


42.0 


B02 


2854 


0.048 


n 


y 


CFHT-Tau5 


G06 


J04325026+2422115 


M7.5 


9.2 


29.8 


G06 


2783 


0.075 


n 


y 


CFHT-BD-Tau 1 


MOl 


J0434 1527+2250309 


M7 


3.1 


19.0 


MOl 


2854 


0.017 


n 


y 


KPN0-Tau9 


B02 


J04355143+2249119 


M8.5 





20.0 


B02 


2642 


0.001 


n 


n 


CFHT-BD-Tau 2 


MOl 


J04361038+2259560 


M7.5 





13.0 


B02 


2783 


0.007 


n 


n 


CFHT-BD-Tau 3 


MOl 


J04363893+2258119 


M7.75 





55.0 


B02 


2747 


0.007 


n 


y 


CFHT-Tau6 


G06 


J04390396+2544264 


M7.25 


0.4 


63.7 


G06 


2818 


0.024 


y 


y 


CFHT-BD-Tau 4 


MOl 


J04394748+2601407 


M7 


2.6 


340.0 


L04,M01 


2854 


0.062 


y 


y 


2MASSJ0455 


L04 


J04552333+3027366 


M6.25 







L04 


2960 


0.015 




y 



2. The TMC brown dwarf sample 

For a typical TMC member, with an age of 3Myr, the stel- 
lar/substellar boundary at 0.075 Mq is between M6 y and M6.5 V 
according to pre-main-sequence tracks from 'Baraf fe et alJ 
( 1998). Therefore, we compiled from the literature objects with 
spectral type later than M6 V to build our bona fide BD sample. ' 

The first optical imaging and spectroscopic surveys of 
BDs in the TMC concentrated on high stellar density regions, 
while the majority of the volu me occupied by the molecu- 
lar c louds was left unexplored (Briceno et al. 1998; Luhmar 
|2000; Bri ceno et al. 2002; Luhman 2 004). Recently, Guieu et al 
(2006) performed a large scale optical imaging survey of TMC 
with the Canada-France-Hawaii telescope with the CFH12k and 
MEGACAM large-scale optical cameras, covering a total area 
of *28 deg^, and encompassing the densest parts of the cloud 
as well as their surroundings. By employing all-sky catalogs, 
iLuhmanI (l2006h considered a 15° x 15° area, large enough to 
encompass all of the TMC. These optical surveys of the TMC, 
conducted since nearly 10 years, yielded the identification of 
42 BDs. 



' In XEST, there are only five objects with a M6V spectral type 
(namely, V410 X-ray 3, MHO 5, MHO 8, KPNO-Tau 14, and CFHT- 
Tau 12), which have luminosities high er than 0.03 Lp, placing them 
above the stellar/substellar boundary llBarafFe et alJl993) . 



2.1. Previous X-ray results on the TMC brown dwarfs 

Two BDs of this sample were previously detected in X-rays by 
the ROSAT/Position Sensitive Proportional Counter (PSPC) - 
which was sensitive only to soft X-rays from 0.1 to 2.4 keV - 
but without being identified as substellar objects^ 

The BD MH04 (B riceno et al. 1998; Luhrnai] 120001) was 
first detected during a 4.0 ks ROSAT/PSPC pointed observa- 
tion of V826Tau in August 1992 with an X-ray count rate of 
4 counts ks ' and proposed as a new weak-lin e T Tauri star 
of Taurus- Auriga (source RXJ043 1.3+ 1800 of Wic hmann eFS] 
(1996). MH04 was also detected during a 7.7 ks ROSAT/PSPC 
pointed observation of the dark cloud Lynds 1551 in February 
1993 with an X-ray count rate of 3. 1 counts ks ' and proposed as 
a new member of the TMC (source L1551X 15 of Carkner et al] 
1996, associated with a "faint, very red star"). Neu hauser et alJ 
(.1999i) reported these two X-ray detections of MHO 4, consid- 
ered at that time as a BD candidate^. 

CFHT-BD-Tau 4, one of the BDs of Martin et al.' ^200)^, 
was detected during a 10.6 ks ROSAT/PSPC pointed observation 
of the Heiles dark cloud 2 NW in March 1993 with an X-ray 
count rate of 4.5 counts ks"' corresponding to a luminosity of 

^ V410 Anon 13, V410X-ray 3, and MHO 5, which were considered 
as BD candidates detected in X-rays by ROSAT { Neuhauser et al. 1 99%, 
have now more reliable spectral types (M5.75V, M6V, and M6V, re- 
spectively; Briceno et al. 2002), which combined with their luminosi- 
ties place them above the stellar/substellar boundary. Therefore these 
very low-mass stars were not included in our sample of BDs. However, 
they were detected in the XEST and will hence be used in ^as TMC 
low-mass stars detected in X-rays. 
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Table 2. Detection list of young BDs in the XEST. Col. (1) numbers correspond to labels in Fig.^ Col. (2), (3) and (4) give the 
BD, satellite and target names, respectively. The naming of XMM-Newton sources in Col. (5) follows the convention of'Giid el et alJ 
(2006a), where the two and three digits code for the field and the source number in this field, respectively. Cols. (6)-(9) give X- 
ray source p ositions, total positi onal uncer t ainties , and distance to the 2MASS position, respectively. References: [FGMSD03] = 
iFavataet alJ 12003) . [BFR03] ^ .Ballv et al.l ( l2003h . There is only pn data for X-ray source #1 as it falls in the gap of the MOS CCDs 
in window mode in this observation. X-ray source #4b is affected by pn gap. For X-ray source #9 only MOS data are available. 



# 

(1) 


BD name 
(2) 


Satellite 
(3) 


Target 
(4) 


X-ray source name 
(5) 


Q'J2000 

(6) 


<5j2000 

(7) 


err. 

(8) 


dist. 

ft 

(9) 


1 


2MASS J0422 


XMM 


TTau 


XEST-0 1-062 


04''22"13!2 


19°34'4a'2 


1.7 


1.5 


2 


KPN0-Tau5 


XMM 


DITau 


XEST- 15-044 


04h29"'45=:8 


26°30'47'.'9 


1.8 


1.9 


3a 


MHO 4 


XMM 


V955 Tau 


XEST-22-021=[FGMSD03] 17 


04''31"24'!2 


18° 00' 21'.' 6 


1.6 


2.1 


3b 


MHO 4 


CXO 


L1551 


[BFR03] 18 


04''31"24':l 


18°00'21'.'4 


0.7 


0.7 


4a 


CFHT-Tau5 


XMM 


GKTau 


XEST-04-003 


04''32"50!3 


24°22'll'.'l 


1.7 


0.7 


4b 


CFHT-Tau 5 


XMM 


V928 Tau 


XEST-03-031 


04''32"'5a;3 


24°22'11'.'4 


1.7 


0.6 


5 


CFHT-BD-Tau 1 


XMM 


CI Tau 


XEST- 17-068 


04''34"U5!3 


22°50'33'.'l 


1.4 


2.5 


6 


CFHT-BD-Tau3 


XMM 


HP Tau 


XEST-08-080 


04''36"38':9 


22°58'13'.'2 


1.9 


1.2 


7 


CFHT-Tau 6 


XMM 


TMC lA 


XEST-05-005 


04''39"'04':l 


25°44'26'.'4 


1.9 


1.9 


8 


CFHT-BD-Tau 4 


CXO 


LI 527 


CXOU 1043947.5-^260140 


04''39"47'':5 


26°01'4a'8 


0.3 


0.0 


9 


2MASS J0455 


XMM 


SUAur 


XEST-26-012 


04''55"23!1 


30°27'38'.'2 


2.0 


3.1 



~ lO^"^ erg s ' (source HCL 2 NW-7a ofBriceno et al.'1999', with 
no opt ical counterpart but a NIR counterpart). Mokler & Stelzer 
(l2002h reported this X-ray source as the counterpart of this BD. 

The higher sensitivity of th e new gene ration of X- 
ray satel lites, XMM-Newton dJansen et alJ ^001) and 
Chandra (^Weisskopf etai] l2002h . allows detection of BDs 
which are from ~ 10 to ~ 100 times fainter in X-rays. 



pie range from a negligible value to about 9 mag with a mean 
and a median of 1.5 and 0.5 mag, respectively. The average low 
extinction of this sample is mainly due to the initial selection 
process for TMC candidate BDs based on optical photometry, 
and therefore biased towards low extinction BDs. 



2.2. X-ray detected brown dwarfs in the XEST 

The XMM-Newton Extended Survey of the TMC (XEST) is a 
large program designed to systematically investigate the high- 
energy properties of young stellar and substellar objects in this 
nearest star-forming region. It is composed of 19 dedicated 
XMM-Newton pointings plus 8 archival XMM-Newton obser- 
vations. The total XMM-Newton exposure tirne is a bout 300 h 
shared among these 27 pointings CGudel et alj|2006a ). We will 
focus here on the X-ray detection of BDs with the three X-ray 
telescopes of XMM-Newton - two European Pho ton Imaging 
Came ras (EPIC) equipped with MOS CCD arrays ( TurrjeretaL 
120011), and the third carrying a pn CCD camera ( Striid er et al. 
[200 1') - simultaneous f/-band observations obtained with the 
XMM- Newton optical/UV monitor are discussed in Grosso et alJ 
S2IM1 . 

Fig. n shows the XEST area and the locations of the 42 
TMC BDs. Fifteen XMM-Newton pointings (of which three 
are archival observations) surveyed 16 BDs. We looked also 
for serendipitous Chandra observations at the Chandra X-ray 
Center's archives (see online Appendix|X| Table lA.fl . We sup- 
plemented the XEST with the ACIS-I observation of L1527 
which allow us to survey and detect one more BD (CFHT-BD- 
Tau4); and with the ACIS-I observation of L1551 which pro- 
vides a second epoch observation of MHO 4. Therefore, we sur- 
veyed in X-rays 17 out of the 42 TMC BDs, which allows us 
to survey the X-ray emission of 40% of the known BDs in the 
TMC. 

Tabled lists the properties of these 17 BDs. Fig. |2] gives 
their location in the HR diagram. The coolest object of the sur- 
veyed BD samp le has a spe ctral type M9.5V, corresponding to 
about 0.02 Mq ( iBaraffe et a l. 2002). Their luminosities range 
from 0.075 Lq to 0.001 Lq. The optical extinctions of this sam- 



M6 M6.5 M7 M7.5 M8 M8.5 M9 M9.5 




3.48 3.46 3.44 3.42 3.40 3.38 
log(WK) 



Fig. 2. HR diagram of the 17 BDs of the TMC surveyed in the 
XEST. The references for the computation of the effective tem- 
peratures and the luminosities are given in Table ^ The pre- 
main-sequence tracks from Baraffe et al. ( 1998) are shown for 
comparison. Continuous lines show mass tracks from 0.1 down 
to 0.02 Mq in steps of 0.01 M©; the dashed line indicates the 
stellar/substellar boundary at 0.075 Mq, which is equivalent to 
spectral type later than M6V for a typical TMC member age of 
3 Myr. A cross indicates typical uncertainties. Two circles have 
been slightly moved in spectral type to avoid overlaps. The 9 
BDs detected in X-rays are marked with 'X' . 
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T he list of X-ray sources in the XEST iGUdel et alJ 
Hop) was completed with X-ray sources of the two archival 
Chandra observations (see online AppendixlAlfor details on the 
data reduction of Chandra data). We cross-correlated the BD po- 
sitions with the positions of all X-ray sources. We detected 9 
BDs in X-rays (see Table |3- The two BDs previously detected 
with ROSAT were redetected with XMM-Newton and Chandra. 
The detection of MHO 4 in LI 551 was previously published by 
[pavata et al. (2003) with XMM-Newton, and Bally et al. (2003) 
with Chandra. Seven BDs are detected here for the first time in 
X-rays. The detection rate in X-rays of these 17 BDs is therefore 
53%. 

We detected mainly the BDs with luminosities greater than 
~0.01 Lq, i.e. the BDs with spectral type earlier or equal than 
M8V (Fig.EJ. 

3. X-ray variability analysis of the TMC brown 
dwarfs 

Obvious variability was observed only from CFHT-BD-Tau 1, a 
nonaccreting M7 V BD, which displayed an X-ray (coronal) flare 
during ~30% of the observation period (Fig.|3- The shape of this 
light curve, a fast rise followed by exponential decay, is typical 
of X-ray flares from young stellar objects (e.g., Imanishi et al. 
l2003t iFavata et alJ[2005l) . We fit the non-zero count rates with 
a quiescent level plus an exponential rise and decay (Fig. |3}- 
We estimate from this light curve fit that the quiescent level 
contributes to 54% of the total counts (quiescent+flare) de- 
tected from CFHT-BD-Tau 1 . The quiescent emission level is 
1.8±0.5 counts ks"', the flare amplitude is 13.3±2.9countsks"', 
i.e. about 8 times the quiescent level. The time scales of the rise 
and decay phase are Trise = 15 ± 5 min and Tdecay = 47 ± 1 1 min, 
respectively. Both values are in the lower range of time scales 
observed in the systematic study of X-ray flares from low-mass 
young stel lar objects in the p Op hiuchi star-forming region with 
Chandra (^Imanis hi et al pool . In this latter star-forming re- 
gion, l^awa et al.l (12005 ) observed with XM M-Newton from 
GY3 10, a young M8.5V BD (Wi lking et al7 T999). an X-ray flai-e 
with a similar amplitude (peaking 6 times above the quiescent 
level) and decay time scale, but with a rise phase about 5 times 
longer. Assuming that the count rate scales linearly with the X- 
ray luminosity (e.g., Mitra-Kraev et al. 2005), we compute, from 
the average X-ray luminosity of CFHT-BD-Tau 1 (see Table 
that the total energy released by the CFHT-BD-Tau 1 flare is 
about 2 X 10^3 ergs. 

The X-ray light curves of the other BDs show no large vari- 
ability (Fig. IB.lt . Two BDs, observed at two differents epochs, 
show variations in X-ray luminosity lower than a factor of 2 
(Fig. IB.2I see below Tables and 0}. We looked for a possible 
modulation of the EPIC light curve of MHO 4 (see F ig. 15721 us- 
ing a Lomb-Scargle periodogram (e.g., Flaccomio et al. 2005), 
but we found no statistically significant periodic signal. The BD 
count rates observed in the XEST are mainly representative of a 
"quiescent" emission, however (micro)flaring activity could be 
present even t hough unresolved in the light curves as observed 
in active stars dCiudel et alJ2()0l . 



4. X-ray spectral properties of the TMC brown 
dwarfs 

4. 1 . X-ray spectra of the brightest X-ray sources 

In total, we have 11 X-ray observations of the 9 BDs de- 
tected in the XEST (Table |2}- We performed spectral fits for 



10 - 



Fig. 3. Light curve of the X-ray flare from CFHT-BD-Tau 1 . The 
bin size is 1000 s. The grey thick line shows the fit of the non- 
zero count rates using a quiescent level plus an exponential rise 
and decay. 



the two Chandra sources, and for these XMM-Newton sources 
with more than ~100 counts (Fig. The computation of 
XMM-Newton source spectra with associated redistribution ma- 
trix files (RMFs) and auxiliary response files (ARFs) is detailed 
in Giidel et al. (2006a), for Chandra observations see online 
Appendix El Spectral fitting was performed with XSPEC (ver- 
sion 11.3: . Arnaud. .1996) usi ng one or tw o optically thin ther- 
mal plasma models (MEKAL; iMewe et d][l995) with 0.3 times 
the solar elemental abundances, combined with an X-ray absorp- 
tion model (WABS; Morrison & McCammon 1983), i.e. the same 
model as use d in COUP to a llow u s to directly compare with the 
result from Preibi sch et alJ ( 12005 b). For XMM-Newton spectra, 
we fitted simultaneously the EPIC pn, MOSl, and MOS2 spec- 
tra. 

Table |3l gives our best fit parameters obtained for 4 BDs 
in 3 XMM-Newton and 2 Chandra observations. The XMM- 
Newton and Chandra spectra of the brightest source, MHO 4, 
are better described with a two-temperature plasma. The two 
plasma temperatures found separately from XMM-Newton and 
Chandra spectra are consistent with each other. For CFHT-Tau 5, 
the hydrogen column density, A^h, determined by one temper- 
ature spectral fitting is 3 (with an acceptable range of 2-4) 
times lower than the value predicted from Ay = 9.2 ± 0.8 mag 
dGuieu et al. 2006) using the A^h Mj ratio o f Vuong et al. (200 j 
combined with the Ay /Aj ratio of lCardeliret al., C1989) for Ry = 
3.1, which leads to A^h - 1-6 x 10^'Avcm"^mag"^^ This may 
point out a non-canonical dust-to-gas ratio on the line of sight 
towards this object. 



4.2. Quantile analysis of the faintest X-ray sources 

For the faintest X-ray sources, we used a quantile analysis, the 
new spectral classification technique for X-ray sources proposed 
by ^ong et al. (2004). This technique uses a quantile diagram 
based on the X-ray colours: x = log((25o/(l - Q50)), and y = 3x 
275/625, where = (E^% - Ei)/(Ef - Ei), with E,,% the energy 
where the net counts are x% of the source net counts in the Ei-Ef 
energy range. We take here Ei = 0.5 keV and Ef = 7.3 keV. 
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^ XEST-Ol-062 XEST-22-021 043134,04+ 180021,0_grpl0.pi 




Fig. 4. X-ray spectra of TMC BDs. From top to bottom and left to right: 2MASS J0422 (pn spectrum), MHO 4 (XMM-Newton and 
Chandra spectra), CFHT-Tau5 {XMM-Newton spectra), CFHT-BD-Tau 1 {XMM-Newton spectra), CFHT-BD-Tau 4 (C/zanc/ra spec- 
tra). Black, red, and green code for XMM-NewtonjEPlC pn, MOSl, and MOS2 spectra, respectively. The continuous lines show our 
best fits obtained with an absorbed single-temperature or two-temperature (MHO 4 spectra) plasma model (see Table|3|l. 



Table 3. Spect ral properties of BDs obtained from spectral fitting. To fit the spectra (Fig. O, we use d a WABS absorptio n model 
JMorrison & M cCammon 1983) combined with one or two MEKAL optically thin thermal plasma model JMewe et alJl995h with 0.3 
times the solar elemental abundances. MEKAL plasma models were computed rather than interpolated from a pre-calculated table. 
We used statistics with standard weighting. Col. (4) gives the net source counts collected by the instrument given in Col. (3). 
Confidence ranges at the 68% level {Ax^ - 1; corresponding to cr = 1 for Gaussian statistics) are given in parentheses. The value 
of reduced and v, the degrees of freedom, are indicated in Col. (10). The emission measures in Cols. (8) and (9) and the X-ray 
luminosity in the 0.5-8 keV energy range corrected for absorption in Col. (11) were computed assuming a distance of 140 pc for the 
TMC. The X-ray fractional luminosity, rf - log{Lx/ L^), is given in the last column. For CFHT-BD-Tau 1, the second line gives an 
estimate of the quiescent X-ray luminosity derived for the light curve fit (see 





BD name 


Instr. 


N 


Nn 


kTi 


kT2 




EM, 




xi (v) 




n 


# 








10-' cm-2 


keV 




10^' 


cm 




10^^ erg s"' 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


(9) 


(10) 


(11) 


(12) 


1 


2MASSJ0422 


pn 


102 


1.7 (0.9-2.8) 


0.7 (0.6-0.8) 






2 




0.83 (07) 


02.0 


-3.5 


3a 


MHO 4 


EPIC 


571 


1.2 (0.7-1.7) 


0.4 (0.3-0.5) 


1.1 (1.0-^ 


1.2) 


5 


6 


0.56 (26) 


08.3 


-3.3 


3b 


MHO 4 


ACIS-I 


271 


(0-0.2) 


0.5 (0.4-0.6) 


1.1 (1.1- 


1.3) 


4 


6 


0.85 (18) 


08.0 


-3.4 


4a 


CFHT-Tau5 


EPIC 


314 


4.9 (3.7-7.7) 


1.5 (1.2-1.9) 






14 




0.82(13) 


11.8 


-3.4 


5 


CFHT-BD-Tau 1 


EPIC 


206 


6.0 (3.9-10.2) 


1.7(1.1-2.6) 






17 




0.59 (06) 


14.3 
07.7 


-2.7 
-2.9 


8 


CFHT-BD-Tau 4 


ACIS-I 


33 


8.6(5.6-13.9) 


0.5 (0.2-0.9) 






27 




0.74 (3) 


24.3 


-3.0 



With the spectral model of ^4.11 we calculated for pn and 
MOS1+MOS2, from the Perl script provided by iHong et all 
(2flQ4)^, a grid of hydrogen column densities and plasma 
temperatures (A^h, kT) in SHERPA'* and the corresponding X- 
ray colours. We adopted on-axis RMF and ARF; and for 
MOS1+MOS2, the MOSl and M0S2 RMFs were averaged 



' The version 1.7 of the ID L and Perl quantile anal- 
ysis softwares of 'Hong et alj ^2003) is available at 
[http : / /hea-www . harvard . e du/ChaKP 1 an e /quant i 1 e r| . 

" SHERPA is a part of the CIAO package. 



and the ARFs added. Fig. |5] shows the resulting quantile dia- 
grams. Grid lines indicate loci of iso-column densities and iso- 
temperatures. 

We noted that the emission lines of the spectral model pro- 
duce folds in the parameter grid, introducing some degeneracy of 
the X-ray colours; i.e. several different values of (A^h, kT) have 
the same £25%, £50%, £75%, and hence the same X-ray colours. 
The degeneracy of the quantile diagram is studied in details in 
online Appendixicl The result of this analysis is the degeneracy 
map plotted in grey levels in Fig. |5] which indicates the num- 




Fig. 5. Quantile diagram for pn (left) and MOS 1 +M0S2 (right). The grid of hydrogen column density and plasma temperature were 
comp uted following the prescriptions of Hong et al. (2004) using in SHERPA an absorption model (WABS; Morrison & McCammo^ 
Il983h multiplied with a single optically thin thermal plasma model (MEKAL; Mewe et al. 1995) with 0.3 times the solar elemental 
abundances. The background map in grey levels indicates the degeneracy level of grid parameters (see online Appendix O. The 
MOS 1 +MOS2 quantile diagram is used only for X-ray sources without available pn data. In the PN quantile diagram, asterisks and 
diamond mark the plasma parameters obtained from spectral fitting; the dashed line indicates the X-ray colour locus corresponding 
to the parameter uncertainties of source #4a (diamond; Tablel^. The spectrum of #3 requires a two-temperature plasma model (see 
TableOll. 



Table 4. Spectral properties of BDs obtained from quantile analysis. Col. (3) gives the instrument name (MOS stands for 
MOS1+MOS2). Col. (4) and (5) indicate the exposure and the net source counts collected by this instrument (i.e. for MOS the 
exposure average and the sum of net source counts). Col. (6) gives the energies below which the net counts are 25%, 50%, and 75% 
of the source net counts in the 0.5-7.3 keV energy range. The resulting position in the quantile diagram (Fig.|5j is given in Col. (7). 
The observed optical extinctions (see Tabled were used to disentangle temperature double solutions in the quantile diagram (#2 
on pn), and to estimate the hydrogen column density (Col. 8) for sources with low constraint in the quantile diagram ( #2, #6 and #7 
on pn, and #9 on MOS), using the relation A^h = 1-6 x lO^'Ay cm^^mag ' jVuong et alJl2003t ICardelh et alJll989ft : these values 
are between brackets in Col. 8. Where negligible optical extinction was measured, we adopted A^h - 0.1 x 10^^ cm"^. When no 
constraint on the plasma temperature was obtained, we adopted 1 keV (value between brackets in Col. 9). The X-ray luminosity in 
the 0.5-8 keV energy range corrected for absorption in Col. (10) was computed assuming a distance of 140 pc for the TMC. The 
X-ray fractional luminosity, rj - log(Lx/Lt), is given in the last column. 





BD name 


Instr. 


Exp. 


N 


Eisvc 1 Eso% 1 £75% 


{x,y) 




kT 


Lx 


1 


# 






ks 




keV 




10^' cm-2 


keV 


10^** ergs"' 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


2 


KPNO-Tau5 


pn 


17.2 


25.3 


0.71/0.82/1.02 


(-1.31,1.21) 


[0.1] 


0.2 


0.8 


-4.1 


4a 


CFHT-Tau5 


pn 


26.1 


118.4 


1.08/1.45/1.89 


(-0.79,1.25) 


6.4 


1.4 


13.9 


-3.3 


4b 


CFHT-Tau5 


MOS 


22.9 


26.6 


1.02/1.27/1.60 


(-0.90,1.42) 


11.1 


0.5 


20.1 


-3.2 


5 


CFHT-BD-Tau 1 


pn 


20.6 


64.7 


1.11/1.36/1.98 


(-0.84,1.24) 


4.5 


1.5 


12.4 


-2.7 


6 


CFHT-BD-Tau 3 


pn 


23.9 


11.6 


0.79/0.98/1.19 


(-1.12,1.27) 


[0.1] 


[1] 


0.6 


-3.7 


7 


CFHT-Tau6 


pn 


14.2 


31.8 


0.92/0.98/1.26 


(-1.13,1.65) 


[0.7] 


0.1 


2.8 


-3.5 


9 


2MASS J0455 


MOS 


123.9 


5.0 


0.75/0.95/1.15 


(-1.15,1.14) 


[0.1] 


[1] 


0.2 


-4.5 



her of different (A^h, kT) values at each position in the quantile 
diagram. 

We calculated for all XMM-Newton sources pn (or 
MOS1-I-MOS2) X-ray colours with background subtraction 
(Table 0} and errors following the prescriptions of Hong et al. 
(12004 . The MOS1+MOS2 quantile diagram is used only for X- 
ray sources without available pn data. We checked for sources 



#1, #4a and #5 that the plasma parameters derived from the 
pn X-ray colours are consistent with the plasma parameters ob- 
tained from spectral fitting marked with asterisks and diamond 
in Fig. |5l For comparison purposes, we indicate in Fig. |5l the 
X-ray colour locus corresponding to the parameter uncertain- 
ties found in the fitting of the XMM-NewtonfEPlC spectra of 
source #4a, the brightest X-ray source (Table The spec- 
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Table 5. Upper limits to the BD X-ray luminosities. The hy- 
drogen column density in unit of 10^' cm"^ given in Col. (3) is 
obtained from the optical extinction using the relation A^e.ai - 
1.6 X Av cm^^mag"' (Yuon^t al. 2003; Cai'delli et al. 1989). 
Exposure times in Col. (4) are for summed EPIC (pn-i-MlH-M2) 
data, in units equivalent for a pn-on-axis observation. Col. (5) 
and (6) give the upper limits at the 95% confidence level for 
net counts in the 0.5-2 keV energy range and the X-ray lumi- 
nosity in the 0.5-8 keV energy range corrected for absorption, 
respectively. The upper limit to the X-ray fractional luminosity, 
J] = log{Lx/Lt), is given in the last column. 



BD name 


XEST 


A'h,21 


Exp. 


N 


Lx.n 


1 




# 


-2 

cm 


ks 


cnts 


erg s"' 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


2MASS J0414 


20 


1.1 


44.7 


115.5 


15.8 


-3.6 


2MASS J0421 


11 


4.8 


38.4 


19.5 


6.1 


-3.3 


KPNO-Tau2 


23+24 


0.6 


22.1 


30.2 


6.0 


-3.6 


KPNO-Tau4 


02 


3.9 


51.8 


25.8 


5.0 


-3.5 


KPN0-Tau6 


14 


0.1 


26.8 


31.1 


6.4 


-3.3 


KPNO-Tau7 


14 


0.1 


25.6 


16.9 


2.5 


-3.8 


KPN0-Tau9 


08+09 


0.1 


51.5 


37.9 


2.8 


-3.1 


CFHT-BD-Tau2 


08 


0.1 


33.9 


25.5 


2.9 


-4.0 



trum of source #3 requires to be fitted by a two-temperature 
plasma model, and has X-ray colours which are outside the grid 
computed with a single-temperature plasma model. Optical ex- 
tinctions were used to disentangle temperature double solutions 
(source #2 on pn; see online Appendix O. and to estimate the 
hydrogen column density for sources with low constrained X- 
ray colours (sources #2, #6 and #7 on pn, and #9 on MOS) using 
the relation A^h_=L6 x lO^'Ay cm"^ mag"' ( Vuon2 et al. 2003; 
ICard elli et al. 1198^ . Where negligible optical extinction was 
measured, we adopted A^h = 0.1 x 10^' cm"^. When there was 
no constraint on the plasma temperature, we adopted kT =1 keV 
(TablegJ. 

Then, in XSPEC we computed from these plasma parameters 
and the source's ARE and RME the X-ray luminosity, corrected 
for absorption in the 0.5-8 keV energy range needed to repro- 
duce the observed net count rate in the 0.5-7.3 keV energy range. 
Eor comparison, the X-ray luminosities obtained by spectral fit- 
ting of sources #5a and #6 and the ones obtained by quantile 
analysis show only a difference of ~0. 1 dex. 

In the following, for BDs observed and detected twice, 
we use the logarithmic average of their X-ray luminosity. Eor 
CFHT-BD-Tau 1 , we use the quiescent X-ray luminosity. 

4.3. Upper limit estimate of tiie X-ray luminosities for the 
undetected brown dwarfs 

We calculated an upper limit to the X-ray luminosity of each 
undetected BD in the 0.5-8 keV band (see Table |5}. Counts 
within 10" of the optical/near infrared position were extracted 
from the summed EPIC soft-band (0.5-2 keV) image. The ex- 
pected number of counts in the absence of emission from the 
BD was determined from the identical region of the correspond- 
ing "reconstructed image", which is composed of background 
and detected sources output by the XMM-Newton Science 
Analysing System (SAS) source detection algorithm, EMLDETECT 
dCiidel etalj l2006 a). An upper limit, at the 95% confidence 
level, to the number of counts in the region from the BD was 
computed using the prescription of Kraft et al. ( 1991), which ac- 
counts for Poissonian fluctuations in the counts from the back- 



ground, and those from the BD^. This results in an upper limit 
to the sum of counts collected by three different detectors. The 
expected number of counts collected by detector /, Ni, from a 
source of flux fx is dependent on the fraction of the total source 
counts that were collected within the extraction region, e;, the ex- 
posure time, fexp.i, and the source energy flux required to produce 
a count rate of 1 counts"', Ki, such that Ni - eiX fexp.i x fx/Ki. 
Therefore, the upper limit to the source energy flux is com- 
puted from the upper limit to the sum of source counts in all 
the detectors, = ZiM, as fx - N/ 2; ei X fexp.iZ-'^i- A model 
source, constructed at the position of the BD from the calibra- 
tion point spread function used by EMLDETECT, was used to cal- 
culate e;, and convolved with exposure images generated by the 
SAS command EEXPMAP to compute the effective on-axis expo- 
sure time at the BD position, fexp.i, on each detector. The on-axis 
response of each detector was generated using the SAS com- 
mand ARFGEN and the appropriate canned response matrix. Ki 
was calculated in XSPEC assuming an isothermal coronal plasma 
of temperature 1 ke V (11.6MK) with metallic ity 0.3 times that 
of the solar corona ('A nders & Grevess3ll989h . and a hydrogen 
column density obtained from the optical extinction using the 
relation A^h ^ 1-6 x lO^'Ay cm^^mag"' (Vuong et al. 20031 
ICardelli et alJll989 ). Upper limits to X-ray luminosities in the 
0.5-8 keV energy range were calculated assuming a distance of 
140pctotheTMC. 

The median of the upper limits of our sample is about 5 x 
10^^ ergs"'. However, the deep archival observation of SUAur 
helped to detect 2MASSJ0455 down to half this luminosity 
level. 

5. X-ray luminosities of the XEST brown dwarfs 
compared to the XEST low-mass stars 

To put the X-ray luminosity properties of the XEST BDs into 
context, we compare t hem to the X-ray lu minosity of low-mass 
stars in the XEST (see lGtidel et al . 2006a). Eor objects observed 
and detected twice, we use the logarithmic average of their X- 
ray luminosity. Eor TMC members that are unresolved multi- 
ple stars in X-rays, we considered the total bolometric luminosi- 
ties. Eig. |6l shows, for the XEST low-mass st ars (L, < 10 Lp , 
which corresponds at 3 Myr to M < 2 Mq, ISiess e t al. 2000|, 
protostars (Class I sources) and BDs, the X-ray luminosities^ 
compared to the bolometric luminosities. Considering only the 
X-ray detected BDs and low-mass (proto)stars, we determined 
with the parametric EM (Expectation-Maximization) algorithm 
implemented in ASURV dFeigelson & Nelsonl ll985) the follow- 
ing linear regression fit: log(Lx/erg s"') = (30.06 + 0.05) + 
(O.98±O.O6)xlog(Lt/L0), with a standard deviation of 0.4 dex in 
log(Lx) for the low-luminosity objects (L,<10 Lq). This relation 
is very similar to the relations found for low-mass stars in other 

5 The BD 2MASSJ0414 is located on the highly structured PSF 
wings of the bright X-ray source V773 Tau, which are not well- 
modeled, and consequently the contribution of V773 Tau in the extrac- 
tion region of this BD is not well estimated. Therefore, we assumed that 
the number of counts observed in the BD extraction region (6" -radius 
here instead of 10" -radius) contains a negligible number of counts from 
the BD, and so the number of observed counts is identical to the number 
expected from background alone. 

* For the low-mass (proto)stars, we use the X-ray luminosities in the 
0.3-10 keV energy band derived from the DEM method in Giidel et ^ 
1 2006a). For a BD plasma with 0.3 times the solar elemental abundances 
and a typical temperature of 1 keV, the difference of X-ray luminosity 
in the 0.3-10 keV energy band and in the 0.5-8 keV energy band is only 
0.06 dex, and hence can be neglected. 
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Fig. 6. X-ray luminosity vs. bolometric luminosity for the TMC 
BDs (diamonds; and arrows for upper limits) and the TMC mem- 
bers detected in the XEST (white and black dots represent low- 
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Fig. 8. X-ray fractional luminosity vs. effective temperature for 
the young BDs of the TMC and single TMC members of the 
XEST. The symbols are as in Fig. |6l The vertical dotted lines 



mass stars and protostars, respectively; Giidel et al. 2006a). The indicates the stellar/substellar boundary. The solid and dashed 



dotted lines indicate from bottom to top an X-ray fractional lu- 
minosity, r] - log(Lx/Lt), of -5, -4, -3, -2. The grey stripe 
shows a linear regression fit (continuous line) and standard de- 
viation (dashed lines) for detected objects with bolometric lumi- 
nosities lower than 10 L©. 



lines show a linear regression fit and standard deviation, respec- 
tively. 



-2.0 



-3.0 - 




M./Mo 

Fig. 7. X-ray fractional luminosity vs. mass for the young BDs 
of the TMC and single TMC members of the XEST. The sym- 
bols are as in Fig.|6l The vertical dotted lines indicates the stel- 
lar/substellar boundary. The solid and dashed lines show a linear 
regression fit and standard deviation, respectively. 



M7 M6 M5 M4 M3 M2 Ml MO 

I 

O 




3.45 3.50 
log(Te„/K) 

Fig. 9. X-ray surface flux vs. effective temperature for the young 
BDs of the TMC and single TMC members of the XEST with 
spectral type MO or later. The symbols are as in Fig.|6l The solid 
and dashed lines show a linear regression fit and standard devia- 
tion, respectively. 
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young clusters jFeigelson &Montmerl3 119991 iPreibisch et aP 
12005 a\ It is consistent with a linear relation between X-ray 
and bolometric luminosity characterized by < log(Lx/Lt) >- 
-3.5 + 0.4 which is valid from the low-mass stars to the sub- 
stellar regime. However, the bulk of the upper limits of unde- 
tected BDs are below this average X-ray fractional luminosity, 
which suggests a lower X-ray fractional luminosity for the BDs 
with about L»<O.O2L0. Taking into account upper limits of un- 
detected BDs, the median of log(Lx/i.) for the XEST BDs is 
-4.0 (see The X-ray fractional luminosity of XEST BDs is 
hence lower than the one of XEST low-mass stars. 

To investigate the relation between X-ray fractional lumi- 
nosity and physical parameters when one moves from low-mass 
stars to the substellar regime, we focus on objects with spec- 
tral type MO or later, which corresponds for an age of 3 Myr 
to masses and luminosities lower than about 1 M© and 0.7 L©, 
respectively. We used the temperature scale given in Table ^ 
to convert M spectral types to effective temperatures. We de- 
termined source masses from effective temperatur e and lumi- 
nosities, using the pre-main-sequence tracks of Baraf fe et alJ 
( 1998). Masses were interpolated linearly along the isochrones. 
For sources located above the 1 Myr isochrone, we prolonged the 
convective tracks vertically. We attributed 0.02 Mq to KPNO- 
Tau4, i.e. the only BD above the 0.02 Mq track. The source 
IRAS S0430 1+261 (spectral type MO), which is located be- 
low the main-sequence was discarded. We restricted our sam- 
ple to stars which are not multiple in order not to introduce ex- 
tra assumptions when distributing the X-ray flux to unresolved 
components. Fig. shows the X-ray fractional luminosity ver- 
sus mass for this sample. Spearman's rank correlation coef- 
ficient computed with ASURV indicates a correlation with the 
probability of the null hypothesis (i.e. no correlation) P{Q) < 
0.06. The EM algorithm yields the shallow linear regression fit: 
log(Lx/iboi) = (-3.5 ± 0.1) + (0.3 + 0.2) x log(M/M0), with a 
standard deviation of 0.5 dex in log(L x/L,). Thi s relation is con- 
sistent with the relation found by Preibisch et al.. (.2005a) for the 
T Tauri stars (M < 2 Mq) of the ONC. This relation imphes that 
the X-ray fractional luminosity decreases by a factor of about 3 
from 1 Mq stars to 0.03 M© BDs. 

Fig. |9] shows the X-ray fractional luminosity versus the ef- 
fective temperature for our sample ranging from ~3840K to 
~2500 K. Spearman's rank correlation coefficient computed with 
ASURV indicates a correlation with the probability of the null hy- 
pothesis (i.e. no correlation) P(0) < 0.11. The EM algorithm 
yields the linear regression fit: log(Lx/iboi) = (-13.1 + 5.4) -H 
(2.7 + 1.5) X log(reff/K), with a standard deviation of 0.5 dex 
in log(Lx/i*). This relation implies that the X-ray fractional lu- 
minosity decreases by a factor of about 3 from hot coronae of 
solar-mass stars to cooler atmospheres of M9V BDs. 

We computed the X-ray surface flux, i.e. the X-ray lumi- 
nosity divided by the surface area of the source, which is com- 
puted from its bolometric luminosity and effective temperature: 
Fx = crT^^Lx/L,, where cr is the Stefan-Boltzmann constant. 
Combined with the decrease of the Lx/L, ratio with the effec- 
tive temperature, this formula implies that the X-ray surface flux 
decreases with the effective temperature with a power-law slope 
greater than 4. Fig. |9] shows the X-ray surface flux versus the 
effective temperature for our sample. The X-ray surface fluxes 
range from ~10'' to ~3 x 10^ ergs"' cm"^. For comparison, the 
X-ray surface flux of the solar corona duri ng the solar cycl e 
varies from -2x10^0-8x10'' erg s ' cm ^ (' Peres et alJ2'00l . 
Spearman's rank correlation coefficient indicates a log-log corre- 
lation with a high confidence level (P(0) - 10 ■*). The EM algo- 
rithm yields: log(Fx/erg s"' cm"^) = (-13.9+4.8) + (5.7+1.4)x 
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Fig. 10. X-ray fractional luminosity of TMC BDs vs. Ha. TMC 
BDs detected in X-rays are marked with 'X' . 

log(reff/K), with a standard deviation of 0.5 dex in log(Fx)- The 
slope of this correlation is well steeper than 4 due to the decline 
of the Lx/Lt ratio with the effective temperature. This correla- 
tion predicts that a M9V BD should have an X-ray surface flux 
in its corona ~ 10 times weaker than in an MO star. However, this 
is still 4 times higher than the X-ray surface flux of the solar 
corona at the solar cycle max imum ( Peres et al. 2004). A similar 
trend was found bv lPreibisc'h et aL(i2005b) for the BDs of ONC 
(see ^for further comparison between the XEST and the COUP 
BD sample). 

6. X-ray fractional luminosity and Ha emission of 
the XEST brown dwarfs 

6.1. Correlation between X-rays and Ha ? 

Studying a BD sample including young BDs (in the ONC, p 
Ophiuchi, and IC348) plus one intermediate BD (TWA5B) 
and one old (LP 944-20) field BD, Tsuboi et al. (2003) found 
a log-log correlation between the X-ray fractional luminosity 
(value independent of the distance assumption), and EW{Ha) 
(logLx/L, = -5.3 + 1.5 X log £'W(Hq')), and noted that BDs 
with £'W(Hci') greater than 100 A were not detected in X-rays. 
We note that including the corresponding (six) X-ray upper lim- 
its, this correlation vanishes. Moreover, this proposed correlation 
is strongly biased by the low X-ray activity of LP 944-20 - de- 
tected only during an X-ray flare (Rutledge et al. 2000) - which 
is likely explained b y the cooling of this old (500 Myr) field BD 
(IStelzer et al.l2006t) . 

EW(lia) measurements are available for all the XEST BDs 
but 2MASS J0455 (Table[T), allowing us to investigate the rela- 
tionship between X-ray fractional luminosity and Ha emission 
from a well-defined sample of young BDs. Fig.llOlshows the X- 
ray fractional luminosity versus £'W(Hci'). We detected in X-rays 
BDs both with low and high £'W(Hci'). We compute Spearman's 
rank correlation coefficient, and find no significant log-log cor- 
relation between the X-ray fractional luminosity and EWiHa). 

6.2. X-ray activity and accretion 

The Ha emission line cannot be used as a tracer of accretion 
in low-mass stars and BDs without a priori knowledge on the 



N. Grosso et al.: X-ray emission from the young brown dwarfs of the Taurus Molecular Cloud 



11 



100 



a 



10 - 



M6 



o 



o 



o 



o 



X o 
o 



o 



M7 



MB 

Spectral Type 



M9 



LO 



Fig. 11. Equivalent width of Ha emission lines vs. spectral type 
for TMC BDs. TMC BDs detected in X-rays ai-e marked with 
'X' . The dashed line shows the saturation limit of chromospheric 
acti vity at los\L(lia)/Lf] - -3.3 (determined in the open clus- 
ters ; iBarrado v Navascues & Martfn.2003.) . 



limit of pure chromospheric Ha emission in these objects. The 
canonical limit of 10 A for EW(lia) was first used to disentangle 
pure chromospheric activity in Weak-line T Tauri Stars (WTTSs) 
and Ha emission line e xcess in Class ical T Tauri stars (CTTSs) 
produced by accretion. lMartinHl998h proposed EW(lla) limits 
depend ing on the spectral tvpe. 'Barr ado v Navascues & Martini 
( l2003h improved this criterion and extended it into the sub- 
stellar domain. This empirical criterion is the saturation at 
log[L(Ha)/L»] = -3.3, based on observations of nonaccreting 
stars in young open clusters, which corresponds physically to 
the maximum amount of energy that can be released in non- 
thermal processes by the chromosphere, i.e. ~5 x 10""* of the 
total emitted energy. Low-mass stars or BDs exceeding this 
limit are accreting. Fig. II llshows ^^(Ha) versus spectral type 
for the BDs of XEST. The dashed line shows the saturation at 
log[L(Ha)/L,] = -3.3 which increases from £'W(Ha)~24A 
at M6 to -148 A at LO JBarrado v Navascues & Mai-tm .200 3). 
MHO 4 and CFHT-BD-Tau3 are just above this Umit, but 
high-resolution spectrum showed no indication of accretion 
(Mohantv et al. 2005). Therefore, we classify these BD as 
nonaccreting. For the other BDs, our classification based on 
^^^(Ha) is in agreement with de tailed high-res olution accr etion 
identification jMuzerolle et alJ20 05b: Mohant v et alJ2005ft . We 
find 6 accreting and 8 nonaccreting BDs. The X-ray detection 
rates of accreting and nonaccreting BDs are 33%+27% (2/6) and 
63%+36% (5/8), respectively^. 

We compute with ASURV the cumulative distributions of 
the X-ray fractional luminosities - including upper limits - 
with the Kaplan-Meier estimator for accreting and nonaccret- 
ing BDs (Fig.ll2>. Two-population statistical methods provided 
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Fig. 12. Cumulative distributions of the X-ray fractional lumi- 
nosities for nonaccreting (continuous line) and accreting (dashed 
line) BDs. In the top panel, diamonds and arrows indicate detec- 
tions and upper limits, respectively of the two samples. Only 
Kaplan-Meier estimator error bars of the nonaccreting sample 
are shown to clarify the plot. Two-population statistical meth- 
ods show that both samples are drawn from the same underlying 
distribution. 
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Fig. 13. Cumulative distributions of the X-ray fractional lumi- 
nosities of XEST (continuous line) and COUP (dashed line) 
BDs. In the top panel, diamonds and arrows indicate detections 
and upper limits, respectively of the two samples. Only Kaplan- 
Meier estimator error bars of the nonaccreting sample are plotted 
to clarify the plot. Two-population statistical methods show that 
XEST BDs are more active in X-rays than COUP BDs. 



' To compute ratio errors we considered that the number of X-ray 
detected BDs and the number of (non)accreting BDs are not exactly 
known, and we combined both Poissonian errors using Gaussian prop- 
agation, i.e. A(a/h) = a/b X ^(Aa/af + (Ab/h)^ with Aa = yfa and 
Ab = yfb. 



by ASURV*^ show that the probability for the null hypothesis that 

Namely the Gehan and Peto-Peto generalized Wilcoxon tests, and 
the Logrank test, whic h are standard methods of univariate survival 
analysis as described bv lFeigelson & Nelsorl <1985l) . 
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both samples are drawn from the same underlying distribution is 
0.75. Therefore, the accreting and nonaccreting BDs have simi- 
lar X-ray fractional luminosities. The median value of the X-ray 
fractional luminosities is -3.9 for the nonaccreting BDs. 

A large difference is found between accreting and nonac- 
creting low-mass stars. In the T MC, the CTTSs ar e 2.2 times 
less luminous than the WTTS (Giideletal. 2006b). A similar 
result was found in the ONC where the CTTSs (selected with 
8542 A Can hne) are 2.1 times less luminous t han the WTTSs 
jFlacco mio e'tal]l2003L iPreibisch et'al]l2005ah . iPreibisch etaP 
(2005a) proposed that in accreting objects magnetic reconnec- 
tion cannot heat the dense plasma in mass-loaded accreting field 
lines to X-ray temperatures. If true, this implies that coronal ac- 
tivity is less affected by the accretion in BDs than in low-mass 
stars. 

The median X-ray fractional luminosity of nonaccreting BDs 
in the XEST is ~4 times lower than the mean saturation value for 
rapidly rotating low- mass (0.22 < M^/Mp < 0.60) field stars 
(log(Lx/-L,) = -3.3: iPizzolato et al.ll2003h . whereas in ONC 
the median X-ray fractional luminosity of nonaccreting low- 
mass s tars is consistent with this saturation level (.Preibisch et al.. 
I2OO5I . 

7. Comparison of brown dwarf X-ray fractional 
luminosities in the XEST and COUP 

IPreibisch et all (l2005bh in COUP studied the X-ray proper- 
ties of y oung BDs spectroscopically identified in the near-IR 

JSlesnick et al 2004) having spectral types between M6V and 

M9V. Eight out of the 33 BDs^ were clearly detected as X-ray 
sources down to a detection limit of log(Lx/erg s"') - 27.3. The 
near-l R selection helped to find extincted BDs ( Slesnick et al. 
l2004l) . Consequently, the apparently low detection rate of BDs 
in COUP is in many cases related to the substantial extinction 
of these BDs. Considering only the ONC BDs with Ay < 5 mag 
reduces the median visual extinction from 5.6 to 2.2 mag, and 
leads to an X-ray detection rate of 40% (6 out of 15 BDs), simi- 
lar to the X-ray detection rate observed for TMC BDs. 

We compute the cumulative distributions of the X-ray frac- 
tional luminosities of the whole sample of XEST BDs, and of 
the COUP BDs with Ay < 5 mag (Fig.[n). Two-population sta- 
tistical methods show that the probability for the null hypothesis 
that both samples are drawn from the same underlying distribu- 
tion is lower than 0.05. The former and the latter samples have 
median X-ray fractional luminosity equal to -4.0 and -4.2, re- 
spectively. Therefore, the XEST BDs are 1.6 times more active 
in X-rays than the COUP BDs. The origin of this difference is 
not yet understood. 

8. Discussion 

We discuss now the origin of the BD X-ray emission in the 
broader context of the X-ray emission of cool stars. Following 
fPreibisch et al. (2005b), we compare the X-ray fractional lumi- 
nosities of TMC BDs with other cool objects. 

Fig.^^shows the X-ray fractional luminosities of TMC ob- 
jects with spectral type equal or later than M5, X-ray detected 
BDs d^eibisch et al. 2005b, and references therein; Stelzer 

' We suppressed from the young BD sample of IPreibisch et alj 
(2005b) COUP 344=HC 722, which was identified as a foreground (old) 
dwarf object in Slesnick et aLl.200.5) 's erratum. This o bject, detected in 
COUP only during an X-ray flare, was also discussed iri lPreibisch et alJ 
42005h) but kept in the reference sample. 



et al. 2006), the evolved BD G1569 Ba b (IStelzei<l2004l). and 
the field BD LP 944-20 dRutledg e et al."'2000'; 'Martin & Bouvl 

2002) . Cool field stars with spectral type M5 or later are 
also indicated: M field stars from Fleming et al. ( 1993); very- 
low mass field stars VB 10 tFleming et al.. .2001). LHS2065 
(Schmitt&Liefke 200^, HC722 (Slesni ck etalJ 12004. 200^ 
Prei bisch et al.il200 5b). All these cool objects show a similar 
level of X-ray fractional luminosities ranging from ~10"^ to 
~10 (note that, compared to Fig.|8] the spectral type range has 
been reduced by a factor of two in Fig.ll4>. 

Based on Baraffe et al. ( 1998) evolutionary tracks, objects of 
spectral type M7 with an age of 1 Gyr are not BDs, but low-mass 
stars twice as massive as a typical TMC BD having an M7 spec- 
tral type and an age of 3 Myr Moreover, such very cool stars 
also have surface gravities about 40 times higher than in a typ- 
ical TMC BD. This shows, as found by Preibisch et al. (2005lJ 
from the COUP BDs, that the X-ray activity of BD coronae is not 
strongly dependent of the BD mass and the BD surface gravity. 
Therefore, this implies that the relation that we have obtained, 
for the young objects of TMC with spectral type M, between 
X-ray fractional luminosity and mass (Fig.^} is in fact the con- 
sequence of the more fundamental relation between X-ray frac- 
tional luminosity and effective temperature (Fig. This latter 
relation agrees with the overall result (field dwarfs and young 
BDs) in Fig. [21 of the 15 sources shown with spectral types 
M8.5V or later, only 4 have any detected quiescent emission; 
the rest are either not detected at all or (in 3 cases) detected only 
during strong flares. By considering a subsample of evolved BDs 
(from ~1 Myr to ~1 Gyr) with a small range of masses (0.05- 
0.07 M0) but with effective temperatures ranging from 3000 K 
down to 1000 K, Stelzer et al. (2006) found evidence for a sim- 
ilar decline of the fractional X-ray luminosity with the effective 
temperature; showing that the atmospheric temperature plays a 
crucial role in determining the level of X-ray activity. 

The growing evidence for a decline of the coronal activ- 
ity with the effective temperature is analogous to the decline 
in chromospheric (Ha) activity seen in field dwarfs at similar 
spectral types, which show a slow decline from mid- to late M 
in Lua/L, and a sharp drop-off around M9V (Mohantv & Bas^ 

2003) . Although a definitive turn-over of the fractional X-ray lu- 
minosity toward L-type BDs is not proven for the present sample 
due to sparse statistics, we note that the fractional luminosity de- 
clines from low-mass stars to M-type BDs, and as a sample, the 
BDs are less efficient X-ray emitters than low-mass stars. We 
thus conclude that while the BD atmospheres observed here are 
mostly warm enough to sustain coronal activity, a trend is seen 
that may indicate its gradu al decline due to the d rop in photo- 
spheric ionization degree dMohantv & BasrJ2003l) . 

9. Summary 

With the XEST we detected 9 (out of 17) young BDs with spec- 
tral type ranging from M6.25 V to M8 V; 7 BDs are detected here 
for the first time in X-rays. This BD sample surveyed in X-rays 
allowed us to investigate the magnetic activity in the substel- 
lar regime. We recovered a well-known relation between X-ray 
and bolometric luminosity for stars (with L,<1OL0) and BDs 
detected in X-rays, which is consistent with a mean X-ray frac- 
tional luminosity < log(Lx/-L.) > = -3.5 + 0.4. For the XEST 
BDs, the median of log(Lx/i'*) (including upper limits) is -4.0. 
The X-ray fractional luminosity of XEST BDs is hence lower 
than the one of XEST low-mass stars. 

A shallow relation is found between X-ray fractional lumi- 
nosity and mass. We show that the X-ray fractional luminosity 
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Fig. 14. X-ray fractional luminosity vs. spectral type for objects of type M5 and later Detections of late M field stars from 
[Fleming et al. ( 1993) are shown as asterisks. The circles show low-mass stars of the TMC detected in X-ra ys (Giidel et al. 2006a|. 
Diamonds and thick arrows show BDs in the TMC. The other X-ray detected BDs (see Preibischet al. 2005b, and references 
therein) are shown by gray filled squares. For very cool objects with strong flares, the values at flare peak are shown by triangles, 
connected by dotted lines to the quiescent emission. Some symbols have been slightly moved in spectral type to avoid overlaps. 



decline by a factor of about 3 from hot coronae of solar-mass 
stars to cooler atmospheres of M9V BDs. Consequently, a rela- 
tion is found between X-ray surface flux and effective temper- 
ature, which implies a decrease of about one magnitude in the 
X-ray surface flux from a MOV star to a M9V BD. 

No significant relation is found between the X-ray fractional 
luminosity and EW(Ha). We used EW{lia) to identify accreting 
and nonaccreting BDs. Accreting and nonaccreting BDs have 
a similar X-ray fractional luminosities. The median X-ray frac- 
tional luminosity of nonaccreting BDs are about 4 times lower 
than the mean saturation value for rapidly rotating low-mass 
field stars. BDs have higher X-ray fractional luminosities in the 
TMC flian in Orion. 

W e confirm, as previously observed in the Chal 
(Utelzer et al. 2004) and in the ONC ( Preibisch et al. 2005b), 
that there is no dramatic change of the magnetic activity at 
the stellar/substellar boundary. Young BDs of spectral type M 
are sufficiently warm to sustain an active corona. The young 
BDs in the TMC, with a median spectral type of M7.5, have on 
average an X-ray surface flux which is 7 times higher than the 
one observed in the solar corona at the solar cycle maximum. 

Deeper X-ray observations of the coolest M-type BDs in the 
TMC are needed to investigate a possible turn-over of the frac- 
tional X-ray luminosity of TMC BDs around spectral type M9 V. 
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Note added in proofs.— 2MASS J04335245+2612548, a TMC BD with spectral type M8.5V (i Luhman e t al."2006') which was 
reported after the acceptance of this paper, is neither detected in X-rays by XMM-Newton/EPlC or in the f/-band by the XMM- 
Newton optical/UV monitor; it is located on the X-ray PSF wings and the {/-band "smoke ring" of ITTau (XEST- 18-030). 
Therefore, adding this newly identified BD to our TMC BD sample wouldn't change our conclusions. 
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Appendix A: Log and data reduction of 
Chandra Observations 

Table lA.ll gives the log of Chandra observations which sur- 
veyed serendipitously the TMC BDs, and used to supplement 
the XEST. 

Data reduction of archival Chandra data were performed 
with CIAO software package version 3.2.1, XSPEC version 11.3, 
the Penn State charge transfer inefficiency (CTI) corrector 
version 1.45, and the acis_extract package version 3.67'". 
We followed Penn State University (PSU) procedures (e.g., 
iGetman et alJ 12005 ^ starting from data reduction from Level 
1 event files provided by the Chandra X-Ray Center (CXC). 
We used the PSU CTI Corrector ( iTownslev et alJl2002r) to cor- 
rect partially the data for CCD CTI caused mainly by radia- 
tion damage at the beginning of the Chandra mission. Source 
detections were performed with the CIAO's task wavdetect. 
2MASS counterparts were used to coiTect small boresight eiTors. 
The acis_extract package was then used to extract source 
photons, estimate local background, construct source and back- 
ground spectra, compute redistribution matrix files (RMFs) and 
auxiliary response files (ARFs), construct light curves and time- 
energy diagrams, and perform automated spectral grouping and 
fitting. 

Appendix B: Quiescent X-ray liglit curves of the 
TIVIC BDs 

Fig. lB.ll shows the X-ray light curves of the brightest TMC BDs, 
which are consistent with quiescent activity. Fig. lB.ll shows X- 
ray light curves of BDs observed at two different epochs. 

Appendix C: Degeneracy of the quantile diagram 

In contrast to the continuum mod els used as illustration of 
quantile diagram bv lHong et al ](f20b4), which connect one sin- 
gle couple of physical parameters to unique X-ray colours, we 
found that the spectra of optically thin plasma (i.e. a continuum 
bremsstrahlung plus emission lines) may have several couple of 
physical parameters producing the same X-ray colours. To use 
the quantile analysis safely, we need to reckon for each physical 
parameter the number of possible values per X-ray colour. 

We divided first the quantile diagram in a regular grid with 
a 0.03 X 0.03 resolution. Then, the uneven grid of X-ray colours 
of each physical parameter was projected to this regular grid us- 
ing spline interpolation. Each physical parameter therefore gen- 
erates a 3-dimensional surface regularly sampled in the X-ray 
colour plane. Finally, we found numerically the number of inter- 
sections between this surface and a line of constant colour, i.e. 
the degeneracy level. 

Fig. IC. II and IC. 21 show the result for the hydrogen column 
density and the plasma temperature grids, respectively, com- 
puted with WABS X MEKAL plasma model with 0.3 times the 
solar elemental abundances and XMM-Newton EPIC pn RMF 
and ARE. The left pannels show the degeneracy maps, which are 
nearly identical when one considers the uneven sampling of the 
model parameter; in Fig.|5lwe plotted conservatively the max- 
imum of these two degeneracy maps. The right pannels show 
for each parameter the corresponding 3-dimensional surface in 
the X-ray colour plane with folds producing multiple solutions. 



Descriptions and codes for CTI correction and acis_extract can 
be found at http://www.astro.psu.edu/users/townsley/cti] 
and http : / /www . astro . psu . edu/xr ay/ doc s/TARA, respectively. 



The quantile diagram is mainly composed of areas where X-ray 
colours correspond to a unique couple of parameters, but areas 
with multiple solutions - from 2 to 5 - are also present. However, 
the area surface decreases with the number of multiple solutions. 

Fig. lC.3l shows an illustration of degenerated solutions \fy- 
1.3. For example if jci < x < X2 there is a unique solution for A^h 
and kT . If X2 < x < xt, there is a double solution: a first solution 
with low A^H and high kT, and a second solution with high A^h 
and low kT . Generally speaking, if there are solutions of A^h 
and kT ordered by increasing values, the solution number n (with 
n ranging from I to A^) is defined as the «*-value of A^h and the 
A^ H- 1 - «*-value of kT . Degeneracy on both parameter can then 
be fixed using extra knowledge on only one parameter. 

We used the knowledge of the optical extinction to de- 
rive the value of the hydrogen column density using the re- 
lation A^H = 1.6 x lO^'Ay cm^^mag ' (Vuongetal. 2001 
Cardelli et al. 1989). The best estimate of the temperature is then 
the temperature solution coiTesponding to the hydrogen column 
density solution which is the closest to this value. 
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Table A.l. Archival Chandra observations which surveyed serendipitously the TMC BDs. Col. (1) gives the n ame of the instrumen t 
used: T (17'xl7' field of view) and 'S' (8.5' x 25.5' field of view) for imaging and spectroscopy ACIS CCD ('Garmire et aH 2003h . 
respectively. Col. (8) give the name of the TMC BD whitin the Chandra field of view. Col. (9) indicates whether we use this archival 
data to supplement the XEST or not (see following notes). Notes: in the ACIS-S observation #3364, the TMC BD KPNO-Tau 2 was 
located 13' off-axis on ACIS-S2, and it is not detected; we obtained here a better constraint on the X-ray luminosity of this BD using 
the sum of the two XMM-Newton exposures, rather than this ACIS-S observation. In the ACIS-S observation #4488, the TMC BD 
2MASS J0421 was located 5.7' off-axis on ACIS-S3, the pipeline detection algorithm found no source at this location. 



ObsID 


ACIS- 


Target 


Nominal pointing 


Start/End obs. 


Exposure 


TMC BD 


Used 
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5]2000 




ks 






(1) 
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(5) 


(6) 
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(9) 
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04''3r32:7 
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2001-07-23T05:10:ll/24T03:52:12 
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MHO 4 
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I 
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2002-12-06T08:30:12/06T14:21:15 


20 


CFHT-BD-Tau4 


yes 


3364 


S 


V410Tau 


04M8"34;6 


28°22'47" 


2002-03-07T06: 16:32 /07T1 1:45:24 


18 


KPNO-Tau 2 


no 


4488 


s 


FSTau 


04'" 22" 00-: 1 


26°58'07" 


2003-1 1-08T12:57:58/08T21:56:34 


30 


2MASSJ0421 


no 
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Fig.B.l. X-ray light curves of BDs in the TMC. From top to bottom and left to right: 2MASS J0422 (only pn data), KPN0-Tau5, 
CFHT-BD-Tau 3, CFHT-Tau 6 (EPIC), CFHT-BD-Tau 4 (ACIS-I), and 2MASS J0455 (only MOS data). On the Chandra light curve, 
black dots indicate the arrival time and energy of individual X-ray photons. These light curves are consistent with quiescent emission. 




Fig. B.2. X-ray light curves of BDs in the TMC observed at two different epochs. Top and bottom row show MHO 4 and CFHT-Tau 5 
light curves, respectively. During the second XMM-Newton observation of CFHT-Tau 5, the source is located on an EPIC pn gap, 
which explains the count rate difference observed between the two epochs. 
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Fig. C.l. Degeneracy map of the absorption column density parameter (A^h) in the quantile diagram. Left: quantile diagram showing 
in the {x, y) space colours the loci of constant A^h values (dotted lines) for a WABS x MEKAL plasma model with 0.3 times the 
solar elemental abundances using XMM-Newton EPIC pn RMF and ARF. Grey levels indicate for each (x, y) value the number of 
corresponding absorption column density values. Regions where a (x,y) value correspond to a unique value of A^h are coloured in 
Ught grey. Right: 3-dimensional shape of the surface generated by A'h from the colour space of the quantile diagram (x,y). Folds 
produce the degeneracy observed in the quantile diagram. 
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^ ^ Degeneracy levels of kT 
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Fig. C.2. Degeneracy map of the plasma temperature parameter (kT) in the quantile diagram. Left: quantile diagram showing in 
the (x, y) space colours the loci of constant kT values (dotted lines) for a WABS x MEKAL plasma model with 0.3 solar elemental 
abundance using XMM-Newton EPIC pn RMF and ARF. Grey levels indicate for each (x, y) value the number of corresponding 
plasma temperature values. Regions where a (x, y) value correspond to a unique value of kT are coloured in light grey. Right: 3D 
shape of the surface generated by kT from the parameter space of the quantile diagram {x, y). Folds produce the degeneracy observed 
in the quantile diagram. 
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Fig. C.3. Illustration of degenerated solutions. Left: horizontal cut of the A^h surface (Fig. lC.ll right) at y - 1.3. Right; horizontal cut 
of the KT surface (Fig. IC. 2 1 right) at the same elevation. The line colour indicates the degeneracy level of the parameter. Numbers 
label the boundary/changing points when one moved on these surfaces along the x direction, and give the correspondance between 
the two parameters. 



